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The modulational stability of the nonlinear Schrodinger (NLS) equation is ex- 
amined in the cases with linear and quadratic external potential. This study is 
motivated by recent experimental studies in the context of matter waves in Bose- 
Einstein condensates (BECs). The linear case can be examined by means of the 
Tappert transformation and can be mapped to the NLS in the appropriate (con- 
stant acceleration) frame. The quadratic case can be examined by using a lens-type 
transformation that converts it into a regular NLS with an additional linear growth 
term. 



1. Introduction 

Intensive studies of Bose-Einstein condensates (BECs) 1 have drawn much 
attention to nonlinear excitations in them. Recent experiments have re- 
vealed the existence of bright 2 and dark 3 ' 4,5 solitons in BECs, as well 
as topological structures, such as vortices 6 and vortex lattices 7 . An in- 
teresting question concerns how such solitary wave structures may arise 
in this novel context of matter waves in BECs. It is well-known that the 
dynamics of the wavefunction in BEC is described (at the mean field level, 
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which is an increasingly accurate description, the closer the system is to 
the zero temperature limit) by the Gross-Pittaevskii (GP) equation which 
is a variant of well-known nonlinear Schrodingcr (NLS) equation 8 . In the 
NLS context, perhaps the most standard mechanism through which soli- 
tons and solitary wave structures appear is through the activation of the 
modulational instability (MI) of plane waves. 

The MI is a general feature of discrete as well as continuum nonlinear 
wave equations. Its demonstrations span a diverse set of disciplines ranging 
from fluid dynamics 9 (where it is usually referred to as the Benjamin- Feir 
instability) and nonlinear optics 10 to plasma physics 11 . One of the early 
contexts in which its significance was appreciated was the linear stability 
analysis of deep water waves. 

The MI has been examined recently in the context of optical lattices in 
BECs both in onc-dimensional and quasi-one dimensional systems, as well 
as in multiple dimensions. In the former case, it has been predicted the- 
oretically 12 ' 13 and verified experimentally 14,15 to lead to destabilization 
of plane waves, and in turn to derealization in momentum space (equiv- 
alent to localization in position space, and the formation of solitary wave 
structures). 

In the present contribution, we discuss the MI conditions for the con- 
tinuous NLS equation (or equivalently for the GP equation) in (1+1) di- 
mensions (1 spatial and 1 temporal) 

iu t + u xx + s|fi| 2 u + V(x)u = 0. (1) 

u in this equation describes the slow envelope complex field dynamics (mod- 
ulating the fast oscillatory dynamics). The subscripts denote partial deriva- 
tives with respect to the index, s G {1,-1} illustrates the focusing (+1) 
or defocusing (—1) nature of the nonlincarity, while V{x) is the external 
potential. 

In section 2.1, we review briefly the results in the absence of the potential 
(e.g., for V(x) = 0). In section 2.2, we examine the case of a linear potential: 

V{x) = -ax, (2) 

which is relevant in experimental situations with gravitational 16 (and po- 
tentially also electrostatic) fields. In section 2.3, we examine the quadratic 
potential of the general form: 

V(x) = -k(t)x 2 , (3) 

which is relevant to contexts in which the (magnetic) trap is strongly con- 
fined in the 2 directions, while it is much shallower in the third one 1 . 
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The prefactor k(t) is typically fixed in current experiments, but adiabatic 
changes in the strength (and in fact even the location of the center) of 
the trap are experimentally feasible, hence we examine the more general 
time-dependent case. Finally, in section 3 we summarize our findings and 
conclude. 

2. Analytical Results 

2.1. No potential 

We start by recalling the results for the modulational stability of the NLS 
(1) without an external potential, i.e. for V(x) = 0: To this end we look 
for perturbed plane wave solutions of the form 

u(x, t) = ((f) + eb) exp[i((qx — tot) + etp(x, t))] (4) 

and analyze the 0(e) terms as 

b(x,t) = &o exp(i/3(x, i)), ip(x, t) — -0 O exp(i/?(x, t)). (5) 

Using f3(x, t) = Qx — Vtt, the dispersion relation connecting the wavenumbcr 
Q and frequency O of the perturbation (see e.g. 8 ) 

{-n + 2qQ) 2 = Q 2 (Q 2 -2s4> 2 ) (6) 

is obtained. 

This implies that the instability region for the NLS in the absence of an 
external potential, appears for perturbation wavenumbers Q 2 < 2s<f> 2 , and 
in particular only for focusing nonlinearities (to which we will restrict our 
study from this point onwards). 

2.2. Linear Potential 

The case of a linear potential is relevant to any context of a costant ex- 
ternal field (gravitational 16 and electric ones being among the prominent 
such examples). In this setting, the NLS is well-known to maintain its inte- 
grate character 17 . Hence, in some sense, we expect that the modulational 
results/conditions will not be modified in this case. 

The simplest way to illustrate that is by means of the "Tappert trans- 
formation" 17 (£ = x + at 2 ) 

u(x, t) — v(£, t) cxp(— iaxt — l-ia 2 t 3 ) (7) 
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(notice however the difference with the expression proposed in 17 ) which 
brings the Eq. (1) back into the form of the regular NLS equation, without 
the external potential, in which case the condition of Eq. (6) applies. 

Hence, the growth terms will now be ~ exp(z(Q£— fit)) with Q = Q r +iv> 
(when the instability condition is satisfied; Q r = 2qQ, cf. Eq. (6)), and 
hence the instability will be developing according to the spatiotemporal 
form: 

u ~ exp (iQ(x + at 2 ) — iCl r t + vt — iaxt — ^-ia 2 t 3 J . (8) 



2.3. Quadratic Potential 

The quadratic potential of Eq. (3) is clearly the most physically relevant 
example of an external potential in the BEC case, given the harmonic 
confinement of the atoms by the experimentally used magnetic traps 1 . In 
particular, to examine the MI related properties in this case, we will use a 
lens-type transformation 8 of the form: 

u(x,t) = r 1 cxp(if(t)x 2 )v((,T) (9) 

where f(t) is a real function of time, £ = x/£(t) and r = r(t). To preserve 
the scaling we choose 8 ' 18 

n = i/e 2 (io) 

To satisfy the resulting equations, we then demand that: 

-f t = Af 2 + k{t) (11) 
-£ T /£ + Aft 2 = 0. (12) 

Taking into account (10) the last equation can be solved: 

e{t)=£{0)exp(Aj /(a)ds^ (13) 

what reduces the problem of finding time dependence of the parameters to 
solution of Eq. (11). 

Upon the above conditions, the equation for v((,t) becomes 

iv T + v« + \v\ 2 v - 2i\v = 0; (14) 

where 



fi 2 = A, 



(15) 
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and generically A is real and depends on time. Thus we retrieve NLS with 
an additional term, which represents either growth (if A > 0) or dissipation 
(if A < 0). 

Eq. (14) allows an explicit solution when v = v(t) (i.e., a spatially 
homogeneous solution). The latter is of the form 

_____ ('.,..,. 2 , ,_ x 



Jo 



exp A(r) + i9 + iAl / cxp(2A(s))ds , (16) 



where A(r) = 4 J Q T X(s)ds, and A a and a are arbitrary real constants. 

A particularly interesting case is that of A constant. Then from the 
system of equations (10)-(12) and (14) it follows that k must have a specific 
form. f,£ and r can then be determined accordingly. In fact, the system 
(10)-(12) and (14) with A constant has as its solution 

k{t) = (t+t*r 2 / w (17) 

f(t) = (t + t*y 1 /8 (18) 
t{t) = 2V2\VTTt* (19) 

r(t) = ln(^)/8A. (20) 

Notice once again, per the assumption of an imaginary phase in the expo- 
nential of Eq. (9), that our considerations are valid only for A £ R. t* in 
the above equations is an arbitrary constant that essentially determines the 
"width" of the trap at time t = according to Eq. (17). 

In this case the modulational condition remains unchanged, but now 
ui satisfies the dispersion relation ui = q 2 — 4> 2 + 2i\, so the growth (if 
A > 0) or dissipation (if A < 0) is inherent in equation (4). Moreover, all 
the terms are modulated by the constant growth (or decay) rate exp(2Ar), 
and the instability (when present) will be developing according to the form 
v r~j exp (i(Q( — O r r) + (v + 2A)t) with O = O r + iv. fl r = 2qQ. 

If k = k(t) is not given by Eq. (17), then A must be time dependent 
(e.g., A = A(i)). Here one cannot directly perform the MI analysis. How- 
ever, still in this case, we have converted the explicit spatial dependence 
into an explicit temporal dependence. An important conclusion that stems 
from this transformation is that the harmonic potential, viewed in the ap- 
propriate frame (of Eq. (9)) can be considered as a form of growth (or 
dissipation, depending on the sign of A for A € R) term. 
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3. Conclusions 

In this brief communication, we have examined the problem of modula- 
tional instabilities of plane waves in the context of Gross-Pittaevskii equa- 
tions with an external (linear or quadratic) potential. Both of these cases 
are directly relevant to current experimental realizations of Bosc Einstein 
condensates. 

It was found that the linear (also integrable) case can be reverted to 
the original NLS setting (wherein the MI conditions are well-known) by 
an appropriate (Tappert) transformation. This transformation was used 
to develop the form of the growth of unstable wavenumbers. Notice, how- 
ever, that these are the wavenumbers of the expansion in a novel spatial 
coordinate which is essentially the spatial variable in a frame with constant 
acceleration. Hence, in this case, in the frame moving with a constant ac- 
celeration (induced by the constant force), the problem resumes its original 
NLS format and the MI conditions and resulting growth can be obtained 
in that frame and then restored in the original frame. 

For the case of the quadratic potential, a lens transformation was used 
to cast the problem in a rescaled space and time frame (in a way very rem- 
iniscent of the scaling in problems related to focusing 8 ' 18 ). In this rescaled 
frame, the external potential can be viewed as a form of external growth. 
For specific forms of temporal dependence of the prefactor of the harmonic 
term (e.g., for appropriate, non-autonomous quadratic potentials), the re- 
sulting prefactor is constant. In such a context once again the MI analysis 
can be carried through completely, producing similar conditions, but now in 
the new dynamically rescaled frame/ variables (which can be appropriately 
re-cast in the original variables). 

It would be interesting to examine if similar considerations can be gen- 
eralized to the case of multiple dimensions. Furthermore, in the current 
scheme of things it seems that the cases of different forms of the potentials 
need to be treated in different ways (which can be understood in terms 
of the different physical effects that they represent). Nonethless, it would 
be very useful, if a general formulation could be developed that could be 
applied independently of the form of V(x), having as special case limits, 
the potential forms presented herein. 

Finally, it would be worth exploring whether the explicitly demonstrated 
as modulationally unstable settings presented herein can be used as a means 
(i.e., initial condition) for directly producing solitary (matter) wave struc- 
tures in BECs in an alternative fashion to the ones currently used in BEC 
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experiments. 
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